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Theory of self induced resonances in asymmetric two-junction interferometer device is presented. 
In real devices it is impossible to have an ideal interferometer free of imperfections. Thus, we 
extended previous theoretical approaches introducing a model which contains several asymmetries: 
Josephson current e, capacitances x and dissipation p presented in an equivalent circuit. Moreover, 
non conventional symmetry of the order parameter in high temperature superconducting quantum 
interference devices forced us to include phase asymmetries. Therefore, the model has been extended 
to the case of 7r-shift interferometers, where a phase shift is present in one of the junctions. 
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INTRODUCTION 



Superconducting quantum interference devices 
(SQUIDs) are the most employed superconductive 
electronic circuits in practical applications.— >2i>^>^>^>^>^>^ 
With the discovery of high-temperature superconduc- 
tors (HTS) also high-temperature SQUIDs have been 
developedjiSiiii'i^'i^iii This class of devices, although 
less sensitive than the most competitive low-temperature 
SQUIDs, have been used in several appHcations, where 
portability and/or positioning as much as high working 
temperatures are needed. Moreover, the demonstration 
of an unconventional symmetry of the order parameter 
in YBaCuO (YBCO),— li^iii opened new horizons 
for using the so-called pi-SQUIDs in superconductive 
electronics. Indeed, tt-SQUIDsAS, can be used to self- 
frustrate quantum bit circuits or to feed RSFQ (rapid 
single flux quantum) devices . As a consequence, 
a full knowledge of properties of HTS SQUIDs is at 
great importance. In particular, the aspects limiting 
their utiHzation in appHcations have to be explored. 
We can consider two effects limiting performance of 
HTS zero-or tt-SQUIDs (zero indicates the conventional 
SQUID where no phase shift has been established along 
the superconducting loop): asymmetries in the junction 
properties and anomalous electrical behavior induced 
by an arbitrary phase shift in one of the two junctions 
forming the interferometer. Asymmetries in conventional 
low-temperature devices have been first examined by 
Tesche and darkest In their paper a complete study 
of the performance in terms of noise characteristics has 
been carried out. The interest on asymmetric SQUIDs 
grew up again after the discovery of HTS. Indeed, the 
parameter spread in HTS SQUIDs is often so large that 
significant asymmetries arise. Hence, it is particularly 
hard to fabricate two identical HTS Josephson junctions, 
even though they are very close to each other on the 
chip. Performance of asymmetric SQUIDs have been 
analyzed by Testa and co-workers .^^'^^ From their 
papers it is evident that higher magnetic sensitivities 
are achieved when asymmetric SQUIDs are used. The 




Figure 1: Theoretical model of an asymmetric superconduct- 
ing quantum interference device and the equivalent circuit 
contains two Josephson junctions with the critical current Ici 
and parallel capacitance d. Each single junction contains 
parallel linear resistance Ri and the interferometer is fed by 
an external source Ic- The self-inductances of the junctions 
are equal Li and ij>i is the phase difference across the ith junc- 
tion. 



asymmetry combined with a damping resistance leads 
to a fiux to voltage transfer coefficient several times 
larger than the one typical for symmetric devices, 
together with a lower magnetic flux noise. The large 
ratio of the flux to voltage transfer coefficient allows 
a direct coupling to an external preamplifier without 
the need of an impedance matching flux transformer or 
additional positive feedback circuitry. This simplifles 
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the read-out electronics, as required in multi-channel 
systems for low-noise measurements. However, the final 
performance of a dc-SQUID is infiuenced by the presence 
of undesired anomalies occurring on the current-voltage 
(IV) characteristics, namely Fiske or resonant stepsj^ 
Such structures originate for the non-Hnear interaction 
between the resonant cavity, represented by the super- 
conducting loop, and an rf current component - the ac 
Josephson current in the junctions. This system may be 
treated with the equivalent electrical resonant circuit, as 
shown in Fig. 1. 

A deep investigation of the properties of resonances in 
asymmetric SQUIDs, also including different phase shift 
in the SQUID loop, is mandatory and can be very use- 
ful for people involved in SQUID design. Self-resonances 
occurring in superconducting interferometers are consid- 
ered to be phenomena reducing performance of high- 
sensitive SQUIDs. Indeed, Zappe and Landman^ first 
investigated experimentally resonances in low-Q Joseph- 
son interferometers. The analysis was taken again by 
Tuckerman and Magerlein;^ who presented a theoreti- 
cal and experimental investigation of resonances in sym- 
metric devices. Successively, Faris and Walsamakis^l 
showed characteristic of resonances in asymmetric two- 
junction interferometers, introducing an important dis- 
tinction between current- and voltage-controlled cases. 
Based on their analysis, Camerlingo et al.— reported 
an experimental work showing the effect of the loop ca- 
pacitance on resonant voltages in asymmetric interfer- 
ometers. Recently, the nature of resonances in SQUIDs 
in which a significant fiux is coupled to the Josephson 
junctions, called spatially distributed junctions (SDJ) dc- 
SQUID, has been analyzed by Chesca.— He showed that 
useful information about the order parameter symme- 
try can be provided by studying directly the magnetic 
field dependences of both the dc Josephson critical cur- 
rent and self-induced resonant modes of dc-SQUIDs made 
of non-conventional superconductors. The further anal- 
ysis of voltage states in current-voltage characteristics 
of symmetric dc tt-SQUIDs, in which the junctions are 
equal and not-distributed circuital elements, has been 
done by Chesca and co-workers.— Moreover, d-wave in- 
duced zero-field resonances in dc tt-SQUIDs have also 
been observed.— 

In our work we present a full investigation of reso- 
nances in asymmetric SQUIDs, also in the presence of 
asymmetries in the junction phases. The outline of the 
paper is the following: In Sec. II we outline the model 
Hamiltonian, and we derive equations for asymmetric dc- 
SQUIDs. In Sec III we present the method and assump- 
tions which have been made. Sec. IV we present our 
results considering special cases and their relevance to 
the other theoretical works. Finally in Sec. V we dis- 
cuss the relevance obtained results to the experimental 
situations. 



II. MODEL 

We start with defining an asymmetric superconducting 
quantum interference device (ASQUID) which consists of 
two Josephson junctions (see Fig. [ij. Each of them has 
a critical current Ici and a parallel capacitance Ci. We 
assume also that single junction contains a parallel Hnear 
resistance Ri and interferometer is fed by an external 
source /c, but the details of the equivalent circuit will be 
specified later. The self-inductances of the junctions in 
ASQUID are equal to Li and L2. We do not consider 
mutual inductances between the junctions. Hamiltonian 
of the ASQUID contains three parts^iiH: 

n = nc + Hj + nM. (i) 

First term on the right side of Eq. |(T|) defines electrostatic 
energy 

nc = lc,v,' + lc,vi, (2) 

where Vi is the voltage across the ith junction. The last 
equation can be transformed to the phase representation 
using the Josephson relation 4> = ^tt/^qV: 

^c^ll^^)\c,<Pl+C,4>l), (3) 

where is a phase difference across the ith junction. 
The second term is the Josephson energy: 

n.! ^ Ejs (1 - cos (t)i) + Ej^2 (1 - cos (t)2) , (4) 

where Ej^i = $o/27r/c,i. To complete the set of equations 
for the interferometer one should take into account that 
loop current II can contribute to the flux. The gauge 
invariant superconducting phase differences between the 
edges of any loop and magnetic flux are directly related 
by the fluxoid quantization relation: 

27r 

02 - 01 = 27rn -I- (pext - -^L+Il, (5) 
<l>o 

where n is an integer and 

L+ = L1+L2, (6) 



Finally, for n — magnetic energy takes the form: 

HM^^L+I.^-^-j . (8) 

At this stage we do not provide an information about 
dissipative environment and external forces which will be 
discussed later. Applying the Euler-Lagrange equation 

dt (d^^C^ -d^^C^O (9) 
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to the Lagrangian 
1 
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2 V27r 

-£^7,1 (1 - COS(/>i) 

we find equations of motion 
— I Ci0i +£;,/iSin0i 

-Z- C202 + i?J,2 Sm (/)2 
ZTT / 



£;j,2(l-COS.^2), (10) 
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Similarly to Tesche^i we introduce the following param- 
eters 
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(1 - x) c, 
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(13) 
(14) 

(15) 



where dimensionless anisotropy quantities x, e and A de- 
scribe the relative deviations of the model parameters 
from the corresponding average values C, Ej and L. We 
can vary the values of the anisotropy parameters within 
the range [0, 1), where zero leads to the isotropic model 
and value one completely rules out presence of one junc- 
tion from the interferometer. Since L+ — L we conclude 
that a difference between inductances does not infiuence 
the dynamics of the model. After renormalization to di- 
mensionless quantities 
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(16) 
(17) 



we finally obtain two coupled non-linear second-order dif- 
ferential equations describing an ASQUID: 



by various excitation sources affects frequencies of the 
oscillating modes of the system. In this paper we as- 
sume that SQUID is current excited by a constant cur- 
rent source (see Fig. [1]). This foundation leads to an 
additional term 7^ in both equations. Origin of the last 
parameter is clear when we consider the equivalent loop 
of a real interferometer— where the center of the induc- 
tance is fed by a gate current source Ig. Using notation 
from Tuckerman's paper and the above information we 
can derive exact form of 7^ : 



In - 2L 
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(20) 



2Ic 2Ic 
where Ic is a circulating current. 

Considering dissipation due to a quasi-particle current 
we add parallel resistances Ri . These dissipative currents 
fiowing through the junctions of the interferometer can 
vary from each other and, as a consequence we have to 
introduce their asymmetry assuming 



il + p)a- 
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(21) 

Different phase shift can be added to each junction 
separately putting (j)i + di in Eq. (fTS]) and Eq. (fT9|) . We 
see that values = and z9i = tt lead to the opposite 
sign of the current which means its opposite direction. 
Finally, we write the equations for ASQUID with phase 
shift in form: 



(1 -t- x) <?i'i + (1 + p) a0i + (1 + e) sin {(f)! + di) 



= 71 + - 

(1 - x) 4>2 + p) a4>2 + (1 - e) sin ( 

{• 

= 72 + - 



62 + 1^2) 

1 - 02) 

13 ■ 



(22) 



(23) 



In order to obtain similar node equations one can also 
use Kirchoff's current law to the specific circuit. We have 
to mention that the noise effects are not present in our 
analysis. Choice of parameters x = e = P = '!^ = stands 
for the fully symmetric case. 



III. METHOD 



(1 - X) '/'i + (1 - e) sin((/.i -f z?i) 



/3 



-,(18) 



(l-x)?i2 + (l + e)sin(</.2 + i?2) = ^^^^.(19) 

Until now we have not considered dissipation effects 
and specific geometry of the circuit. First, we have to 
decide, what modes of operation we think about: cur- 
rent controlled (CC) or voltage controlled (VC)? This is 
a crucial point simply because a choice we make is going 
to affect our system. For the VC case where SQUID is 
excited by a voltage source Vg we have to add terms pro- 
portional to Vst to the equations. The difference caused 



We shall analyze two coupled differential equations 
(|22l) and (|23|) for the case (3 < I that coupling between 
the two junctions of the interferometer is strong and, 
hence the last terms of the right hand in Eqs. (|22|) and 
(|23|) play important role since they contain expressions 
proportional to ±(3~^ (02 — <f>i)- Let us introduce new 
variables 
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where (f)- represents the flux number (this parameter dis- 
tinguishes interferometer from a point junction), and 0+ 
is the average phase difference of the junctions. For the 
equivalent circuit of a real interferometer 7_ can be rec- 
ognized as a control current Ic and 7+ as a bias current 
Ig . Parameters i?± are relative changes of the phase shifts 
present in each junction. In terms of the above we write 
equations (|22|) and (|23|) in form 



-i?+)cos(0_ -7+ 
>_ +i?_)cos(0+ 

= 0, (27) 



+X4>- + ap<t>- + e sin 



P 

+ apnuj + eJ„ {S) cos {9 - 1)+) sin ((^^ + ^36) 



-xSlo^ = J:,7((5)cos(6I~?9+)cos((/.c + i?-) 

+ eJ-{d)sm{e-'d+)sm{(t)c + 'd^), (37) 

apSu) = -J+(5)sin(6'-??+)cos((/)c + ^^-) 

+ eJ+(5)cos(0-7?+)sin(0, + i?_), (38) 



(/)- + a(j)- + sin ((/)_ + i9_) cos (0+ + - 7- + —0- 

+ ap</>+ + e sin (0+ -I- i?+ ) cos ((/>_+ ) 

= 0(28) 

In the following we have to assume a form of the solution. 
The voltage variations appearing in ASQUID come from 
the interaction between the junction current and circuit's 
elements. We assume voltage sinusoidal variations with 
dc component V, ac amplitude v, frequency to and phase 
ip: 



V (t) = V + v cos {ujt + (p) , 



(29) 



where other harmonics are filtered out. Using the Joseph- 
son relations and integrating out we get for zth junction: 



(j>i (t) — ipQ^i + tot ±6 sin {ut + ip) 



(30) 



where S = ^. The flux number (f)- and the average phase 
difference can be written: 



b- ~ 4>c — 5 sinijjt, 
6+ ~ ncut — 6. 



(31) 
(32) 



where (j)c is the average value of the internal phase 0_ . 
In order to account the difference between odd and even 
behavior of the ASQUID interferometer we deflne: 



t>- — (f>c ~ 6 sin Lut ~ k- 

TT 

^ 2 



(33) 
(34) 



where k is equal (1) for even (odd) number of reso- 
nances. 



IV. RESULTS 



Substituting expressions ([331) and ([341) into equations 
([27l) and l[28|) and extracting by calculating average over 
time the dc, smut and coswi Fourier components we get: 

anu} = 7+ — J„ ((5) cos (6* — z9_|_) sin (0c + ??-) 

-t- eJniS)sm{e-i}+)cos{<j)c + ^-), (35) 



5 = J- {S)sm{0~d+)sm{(f>c + d-) 

- (S) cos {9 ~ 1?+) cos {4>c + i3-)(39) 



aSuj = J+ iS)cos{6 ~d+)sm{<j)c + ^-) 

- eJ+ {6)sm{e -^+)cos{4:c + ^-) , (40) 



where 



Jt (S) = Jn-l (S) ± Jn+1 (S) 



(41) 



and J„ (S) is the Bessel function of the flrst kind.= Using 
Eq. ([35]) , Eq. ([40)) and Bessel function identity, we obtain 



(42) 



anuj = 7+ 1 , ^ aou) 

J,t [5) 



1+ 



2n 



We define normaHzed excess current due to the resonance 



2n 



(43) 



Above equations can be rewritten using the dimensionless 
damping parameter T = (aw^) , where Ur is the reso- 
nant frequency. Gamma was introduced by Werthamer— 
and described the strength of the coupHng of the current 
to the resonance in case of the junction coupled to cavity. 
Several authors used it as a damping parameter i^^^^^ We 
can combine equations ([36]) and ([M)) : 



(44) 



2 , apS'^uj 

7- = -5<Pc + apnuj + — , 

p 2n 

and using relation for excess current we get 



7- = ^s'^c + P7+- 



(45) 



We see that formula ([43|l derived for the excess current is 
universal in such sense that it holds even for asymmetric 
SQUID. This expression is also true in the presence of 
any changes of the phase shift in one of the junctions of 
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Figure 2: Current voltage characteristics with Josephson cur- 
rent anisotropy e, first resonance, F — 20, /3 = 0.1. Black 
color of the curves used in this and next plots indicates the 
symmetric SQUID x ~ ^ ~ P ~ 0. 



Figure 3: The normalized resonant current lexc versus damp- 
ing parameter F for different values of the Josephson current 
anisotropy parameters e (red curves), nth resonance. 



to form 



We can add squares of expressions 
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(46) 



where Cj = oj / ujr = V/Vr is the normaHzed voltage. From 
above and Eq. (|43l) we can derive normaHzed excess cur- 
rent dependence on voltage for given anisotropic param- 
eters. However analysis is complex and is better to sim- 
plify our model considering special cases which could give 
us more insight into structure of resonances in ASQUID. 



A. Special cases 

For a general choice of parameters equations (|35l) - l(40|) 
are coupled and must be solved numerically. However 
considerations of special cases can provide more insights 
into general solution of the problem. 



(47) 



From the above expression coupled with Eq. I|43p we can 
derive the normalized current dependence on normalized 
voltage plots with e asymmetry (see Fig. [21). Also the 
normalized resonant current versus damping parameter 
for several resonances can be obtained (see Fig. [3]). 



2. Asymmetry of the capacitances (x 7^ 0) and resistances 
(p 7^ 0) with phase shift {'d± / 0) 

Let us consider case when e = which means that 
asymmetry of the Josephson current is not present. In 
this case equations l(35 l) - ()40| are reduced to 

anuj = j+- Jn{S) cos {6- d+) sin {<f)c + ■&-), (48) 
9 

1?+)C0S (0C+??-) 



-;(/'c - J,i (S) sin (6* 

P 

+apnu!, 



xSuJ^ = J„ ((5)cos( 



(49) 

i)+)cos{cbc + ^-), (50) 



1. Asymmetry of the Josephson current (e 7^ 0) 

In that case we assume that only Josephson current 
asymmetry is present. Then the Eq. (|46l) can be reduced 



ap6uj = -J+ ((5) sin (61 cos (0c + (51) 



J- ((5) sin (6* - i?+) sin ((/.^ + (52) 
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Figure 4: Dependence of the normalized resonant current lexc 
with relative phase shift — 0.5 (red), 0.25 (dot), (black) 
and the ratio T/F with no phase shift &- = versus magnetic 
field (be- 



a6uj = J+ ((5)cos(6l-i?+)sin(<?!)c + ??-). (53) 

First we will analyze low-F case in order to compare 
our results with the original calculations presented in 
literature . ^^1^^ 



Low-V devices with symmetric values of the Josephson 
current (e = 0) 

At the resonance frequency oj = nuJr Eq. l(52l) is satis- 
fied when 9 = 'd+. This condition rules out equations 
with terms proportional to sin(6' — and therefore 
there is no trace of the asymmetries of the Josephson 
current e and dissipation p. For small gamma F devices 
{S) = for n > 1 and, hence only the first resonance 
exists. We can derive the following equations 



S_ 
T 



= cos ( 



- = sm( 



(54) 
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where T = {x^^) ^ is dimensionless parameter. Rear- 
ranging the last equation and putting into expression for 
excess current we get: 



lexc = r sin 



^_) 



which is general result for different SQUIDs. 
calculate other relations: 



(56) 
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Figure 5: The normalized resonant current ^cxc versus 
damping parameter F for different values of the dissipation 
anisotropy parameters p (red curves), nth resonance. 
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Figure 6: Current voltage (/cxc — V/Vr) characteristics with 
dissipation anisotropy p, first resonance (n = 1), F = 20, /3 = 
0.1. 



lexc = — COf 
T 
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which are plotted in Fig. [H We see that the results 
obtained previously by other authora^^*^^ are presented 
in framework of our rather general calculations and can 
be derived as special cases. 



Analysis for not small F 



When F is not small we cannot simplify equations us- 
ing condition under which Bessel functions can be ap- 
proximated by zero except the case of the first resonance. 
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Figure 7: The normalized resonant current Icxc versus damp- 
ing parameter F, different capacitance anisotropy parameter 
X, first resonance {n — 1). 



Putting e = in Eq. (|46|) we obtain: 
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Jn (S) 







+ 





2 

+ 


apSuj 


.Jn{S)_ 




[jni6)\ 



= 1. 



(59) 



From the above equation and expression l|43p for the ex- 
cess current we can derive the normaHzed current voltage 
characteristics for ASQUID. The second and the fourth 
terms of above equation can be combined. We see that 
influence of the anisotropy of the dissipative current 



[1+p] 



aSuj 

Jrt (6) 



(60) 



manifests by the decreasing of the maximum value of the 
resonant current, for given nth resonance mode, when we 
increase the anisotropy parameter p (see Fig. [5] and Fig. 
[6]) . We observe a shift of the maximum value of lexc to- 
ward higher values of the damping parameter T. Analysis 
of the influence of the anisotropy of the capacitances can 
be done in the same manner. We can again merge first 
and third terms of the Eq. l(59l) . Contrary to previous 
simple case present one is more complex merely because 
we have taken into account element proportional to a)^ 
which produces minor changes (see Fig. [7] and Fig. [HI). 
Now even small deviations of the anisotropy parameter x 
from equilibrium have a major impact on equations and 
in consequence on behavior of the ASQUID. For small 
values of X, at fixed value of the damping parameter T 
there are two possible solutions even for the first reso- 
nance. In symmetric SQUIDs this situation was present 



Figure 8: The normalized resonant current /cxc versus damp- 
ing parameter F, for different values of the capacitance and 
dissipation anisotropy parameters x {red) = p (blue), first res- 
onance (n = 1). 



for higher resonances n > 3. Explanation of the lat- 
ter comes from the fact that the resonant circuit oscil- 
lates at a frequency oiojr, while Josephson current in the 
junctions oscillates at ncjr- In ASQUID we have three 
natural frequencies^! 2 = (L+C1.2) and related 



2 which can be excited by the ac Josephson 
effect and converted through nonlinear interactions be- 
tween junction and resonant circuit into dc current steps. 
Therefore introducing capacitance anisotropy we are able 
to create higher modes multivalued behavior of the excess 
current even for the first resonance. 



3. Symmetric case (x = ^ = P — 0) with phase shift 
/ 0) 

This case corresponds with a situation where different 
phase shift is present in the junctions of the interferom- 
eter and analysis is similar to one carried by Chesca.— 
The equations take form: 



anu} = 7+ — J„ ((5) cos (6* — i?+) sin (0c + I?-) , (61) 



7- = -0c- Jn(^)sin(e'-?9+)cos((/.c + t?_), (62) 



J- {5) sin (6* - 7?+) sin ((/-c + d-) (63) 



a5uj = J+((5)cos(e'-i9+)sin(0c + 'i?-)- (64) 
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Figure 9: The normalized resonant current -^exc versus nor- 
malized voltage V/Vr, with different values of the anisotropy 
parameters, capacitance x, the Josephson current e, and dis- 
sipation p for second resonance {n — 2). Black curves refer to 
absence of anisotropy parameters. 



We do not expect any changes in excess current - voltage 
characteristics. Rather, as it was pointed out by Chesca 
the difference between SQUIDs with various phase shifts 
can be visible only in magnetic field. In order to calcu- 
late excess current dependence on magnetic field we add 
squares of the equations (|63l) and ((64l) . The resonant 
current is maximized when 9 = d+ and we can write the 
solution in parametric form: 



[/ea;c;sin {iKpe +"&-)] 



2Tn' TJ+{5) 



(65) 



where (5 is a dummy variable. Changing value of the 
parameter d- from to — 7r/2 we have — and — tt 
SQUID respectively. The shape of the surface describe 
current magnetic field dependence (see Fig. [lO]) remains 
unchanged but is translated by a vector [0; — along 
(be axis. 




Figure 10: Excess current lexc versus magnetic field 0e (first 
resonance, n = 1) characteristic for different values of the 
damping parameter F for — tt interferometer. 



V. DISCUSSION 

The resonances in SQUIDs are investigated theoreti- 
cally with several asymmetries: Josephson current e, dis- 
sipation p and capacitance x- In real devices it is im- 
possible to have an ideal interferometer free of imperfec- 
tions. In practice various deviations of the interferometer 
parameters from average values can occur together and 
mutually conceal each other. At this stage we have to 
separate discussion related to low- and high-Tc SQUIDs. 
In the former case, experimentally, we are able to control 
asymmetry of dissipative parameter p adding a parallel 
resistor to the junction but it is difficult to change the 
Josephson current independently from the capacitance. 
To produce the asymmetry of the Josephson current in 
the interferometer we can change the area of the junc- 
tion A or thickness of the barrier d. Parallel-plate capac- 
itor with area A of the plates and space d between them 
has the capacitance equal C — CrtoA/d for A » c?^, 
where er is the relative dielectric constant of the inter- 
layer dielectric and eo is the vacuum electric constant. 
On the other hand the critical current can be written 
as Ic = jcA where jc is the critical current density. 
These two simple relations imply that varying area AA 
of the junction in the interferometer we change both ca- 
pacitance and critical current proportionally Ale ~ AC 
at the same time. When no further resistor is added to 
the junctions not only capacitance and Josephson cur- 
rent are related. From Ambegeokar-BaratofS^ formula 
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aS5 1.0 l.OS v/V^ O.iS 1.0 l.OS v/V^ 

Figure 11: Current voltage characteristics (/cxc — V/Vr) for 
several asymmetric configurations of the SQUIDs related to 
the Ambegaokar-Baratoff formula (A/c ~ AC ^ AR'^) for 
changes of the junction area AA, first resonance n — 1, T = 
20, /3 = 0.1. 

we know that the product IcRn, where Rn is the re- 
sistance in normal state, has an invariant value which 
depends only on the material in fixed temperature. Thus 
changing the value of the Josephson current we alter the 
resistance of the junction. Recapitulating these rather 
simple considerations we can introduce asymmetry in 
the Josephson current changing the area of the junc- 
tion {Ale ~ AC ~ AR~^). Setting parallel resistor 
we can control value of the resistance and vary dissipa- 
tive parameter independently from the current asymme- 
try. We can also imagine junctions with different thick- 
nesses of the barrier but technologically this case is dif- 
ficult to achieve thus we do not consider it. In experi- 
ments with ASQUID both technically reached asymmet- 
ric cases do not differ very much because of the capaci- 
tance anisotropy. As we see from Fig. [11] the biggest im- 
pact on the maximum value of the resonant current has 
the anisotropy of the capacitance. Even small changes of 
X can decrease excess current almost to zero. 

The situation changes completely when high-Tc 
SQUIDs are considered. On one hand, the probability 
to find junction parameter asymmetries is particularly 
high, because high-Tc junctions are intrinsically affected 
by defects, as for instance faceting and/or oxygen vacan- 
cies inside the barrier. Moreover, up to now, the charge 
transport process is not completely understood, although 
various hypothesis have been proposed,i2L2£iii and other 
recent experiments are still in progress i^^i^° In particular, 
the simple rule IcRn — const vaHd for low-Tc SQUIDs 
does not apply in the case of high-Tc interferometers 



typically used in applications, based on the symmetric 
bicrystal c-axis [001] devices, and changing one single pa- 
rameter is now possible. In such interferometers, IcRn 
is proportional to the critical current density Jc at low 
values and stays roughly constant at high- Jc values 
Moreover, HTS junctions are intrinsically shunted and 
SQUIDs are fabricated with no additional shunt resistor. 
As a consequence, the way to fabricate HTS SQUIDs 
with symmetric junctions is probably to reduce junctions' 
widths, limiting the effect of the interface defects. In all 
other cases, asymmetries will be very probable and our 
analysis could be relevant to understand the presence of 
resonance steps. 

Different approach is necessary in the case of asym- 
metric [001] or [100] HTS bicrystal junctions, where the 
relation IcRn seems to be similar to the one of low- 
Tc systems^ and the necessity to account for effects 
of a non-conventional symmetry of the order parameter 
forces to include also the phase asymmetries in study- 
ing dynamical states in HTS interferometers. Finally, 
also the inclusion of the second harmonic term in the 
Josephson current in order to account for experimental 
resulta^^^"^ is mandatory. This will be the argument of 
a separate paper, and the possibility to deal with one sin- 
gle asymmetric parameter is now eventual. Moreover, a 
non conventional symmetry of the order parameter forces 
to include also phase asymmetries in studying dynamic 
states in high-Tc interferometers. In this frame the cal- 
culations derived in the present paper allow to investigate 
SQUID dynamics in both low- and high-Tc asymmetric 
devices. 



VI. SUMMARY 

In this paper we have presented a detailed theoretical 
study of the resonances in the asymmetric superconduct- 
ing quantum interference device. Analytical approach 
revealed the nature of the resonances in the presence of 
several asymmetries: Josephson current e, capacitances 
X and dissipation p. Also we were able to derive magnetic 
field dependence of the excess current in presence of the 
magnetic field and phase shift. Our calculations imply 
that deviations of the capacitances from the average value 
in SQUID have profound impact on physics of the system. 
We have found that our theory can be useful to deter- 
mine asymmetry parameters present in lightly damped 
ASQUIDs. Especially for SQUIDs produced from HTS 
materials where deviations from average values are prac- 
tically inevitable our considerations are very helpful. 
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